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In this Supplementary we review the superconducting symmetry in nanowires with strong 
spin-orbit coupling as a comparison to the discussion on the 3D TI in the main text. A 
single spin branch is now only achieved after applying a Zeeman magnetic held. We review 
in which regime a dominant p-wave correlation is present and under which conditions a 
dominant p-wave correlation can lead to a zero-energy Majorana bound state. 


I. PAIRING WAVE FUNCTION AND MAJORANA-MODES IN NANOWIRES 


The chemical potential is usually close to the conduction band so that p 3> A. If we also 
consider low energy excitations we obtain for the nanowire {k^ = 0, = A) 
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(1) 

( 2 ) 


The energy dispersion relation of the s-wave proximized nanowire can be obtained by 
diagonalizing the corresponding Hamiltonian as described above, and is found to be 

E = ±\JEt± 2 a/A^, (3) 

where Et = A'^ + el + A"/c" -|- M^ and Eg = elX^k"^ + A^M^ -|- e^M". In the limit of p 3> A 
this can in good approximation be written as 

E = ±ep± VX^k^ + M2, 

to calculate the pairing symmetry in three different regimes M a/A^ -|- p^, M = 
a/A" -I- p" and M 3> a/A^ -|- p^. 
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(a) (b) 



(c) 



FIG. 1. (a) Energy dispersion relation for a nanowire with M <C \/The solid line 

corresponds to the relation Et — 2t/E 7 and the dashed line to Et + (b) M = ^/S? + lE 

the gap at /c = 0 is closing. Because of the finite magnetic field there is another gap indicated by 
A 3 , (c) In the regime M + fj? a dominant topological non-trivial p-wave gap exist at Ai. 

The graphs are plotted for A in the order of p for clarity. 


A. Regime 


In this case we have four superconducting gaps as indicated in Fig. 1(a). Ai and A 2 
indicate the position for positive wave vector momentum. At Ai (and for the corresponding 
gap at negative k) it follows from the dispersion relation that fi = — X\kF\ and for A 2 

/i = - 1 ^ + Xlkpl- Using the relations for the anomolous Green function, Eq. (1), we get for 
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Ai 

~ -2AikkfX^/Znw, 

~ -‘^^X^k'^/Znw, 

F^i ~ 2AX^k^/Znw, 

Fjj ~ 2AikkfX‘^/Znw 

Expressions such as are set to zero in the above equations as we consider low 

energy excitations. For A 2 we obtain 

~ 2AikkfX‘^/Znw, 

F^j ~ -2AX^k‘^/Znw, 

Fji ~ 2A}?k^jZnw, 

Fjj ~ —2AikkfX‘^/Znw 

We see also here that in the case of time-reversal symmetry there is an equal admixture of 
s and p-wave correlations. 

B. Regime M = \/A^~T7? 

This case is shown in Fig. 1(b), where the superconducting gap at /c = 0 is closed. 
At A 2 the anomolous Green function parts are zero because k = 0. Because we have 
turned on a magnetic held, the bands corresponding to different helicity are not orthogonal 
anymore. Therefore, there exist an interaction between the two in this regime and another 
superconducting gap is opening above /x (A 3 in the Fig. 1(b)). [1] For A 3 the relation jj, = 
holds so that 


Fjj ~ 2AlkXfxfZfi^^ 

Fjj ~-A (k^X^-lJ.^) /Zr,^, 

Fjj ^ a jZn^, 

Fjj ~ 2AikXll/Znw 

For a typical [2] spin orbit strength of A = 0.2 eV-A and m = O.OlSrrie where me is the 
free electron mass, it follows that the magnitude of the spin-singlet order parameters in 
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this regime are larger than the p-wave components. However, since the energy level of A 3 
increases with M, i.e. ~ V for low energy excitations this snpercondncting 
contribntion can be neglected. 

^ 1 ’ ~^2rn~l^~ a/= 0 SO that we get the following relations for the anamolons 
Green fnnction part: 


-^11 ~ —2A.ikX X- I 

Fu ~ —2^\^k‘^/Znw, 

F|i ~ 2/^\^k‘^/Znw, 

Fjj ~ 2Mk\ (/r - - /i2) 

Here, we see that Fjj becomes gradnally smaller and F|i becomes larger than the s-wave 
pairing parts, assnming M to be positive. We obtain a dominant p-wave component for the 
spin triplet component where the spins are parallel to the magnetic held. 


C. Regime M > VA2 + p?, “Majorana” regime 

The sitnation is sketched in Fig. 1(c). The gap at /c = 0 reopens. The pairing wave 
fnnction at Ai satishes 


Fii ~ -2Aik\ [m + ^\^k) + M2) /Zn^, 

F|j ~ -2/X\^k‘^/Znw, 

F|i ~ 2lS.}?k^ jZnw^ 

Fjj ~ 2^.ik\ [m - ^\^k) + M2j /Znv,. 

We see that for increasing M the p-wave part becomes more prononnced and that Fjj becomes 
larger than Fjj. Fjj actnally goes to zero for increasing positive held. For negative held Fjj 
will be dominant as one would expect for the opposite held direction. Around k = 0 there 
is no interaction anymore between the hole and electron branches and, as noted before, the 
degeneracy in the p-wave order parameters disappears in this regime. 
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At A 3 we have 


F 11 ~ —2Aik\M/Zmi,, 

Fij ~ 2A (M^ - /Zn^, 

Fji ~ -2A (m2 - k^\^) jZ^^, 

Fjj ~ 2Aik\M/Znw 

Hence, the s-wave pairing becomes dominant for increasing magnetic held at A 3 , but as 
discussed before it is safe to ignore it in the low excitation regime. 

Since the gap at /c = 0 closed for M = a/A^ + and has now reopened again we also 
have a proper localized single Majorana mode at the ends of the nanowire. Just as the 
surface states of a topological insulator are localized to a surface, the Majorana mode in 
a nanowire is localized at the ends of the nanowire due to the vacuum gap at one side 
and the superconducting gap with a reversed band order on the other side. The reversed 
band order ensures that the bands are crossing at the ends of the nanowire resulting in a 
localized end mode similar to the topological surface and edge states in topological insulators. 

Note, that all of the expressions in this section can be generalized to two dimensions. 
This could be relevant for InAs thin dims for example. The scalar k has to be replaced by 
a k-vector with an amplitude |k| and phase e*® analogous to the 3D topological insulator 
discussed below. The e*® factor denotes the chirality of one branch and the chirality of 
the other branch present at the Fermi level. Apart from an additional e*® and factor in 
the expressions of the order parameter, the given conclusions would still be valid. 


[1] P. San-Jose, J. Cayao, E. Prada, R. Aguado, New Journal of Physics 15, 075019 (2013). 

[2] V.Mourik, K.Zuo, S.M.Frolov, S.R. Plissard, E.P.A.M.Bakkers, L.P.Kouwenhoven, Science 12, 
1222360 (2012). 
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Abstract 

To guide experimental work on the search for Majorana zero-energy modes, we calculate the 
superconducting pairing symmetry of a three-dimensional topological insulator in combination 
with an s-wave superconductor. We show how the pairing symmetry changes across different 
topological regimes. We demonstrate that a dominant p-wave pairing relation is not sufficient to 
realize a Majorana zero-energy mode useful for quantum computation. Our main result is the 
relation between odd-frequency pairing and Majorana zero energy modes by using Green functions 
techniques in three-dimensional topological insulators in the so-called Majorana regime. We discuss 
thereafter how the pairing relations in the different regimes can be observed in the tunneling 
conductance of an s-wave proximized three-dimensional topological insulator. We discuss the 
necessity to incorporate a ferromagnetic insulator to localize the zero-energy bound state to the 
interface as a Majorana mode. 

PACS numbers: 74.78.-w, 74.78.Fk, 74.45.-be, 03.65.Vf 

Keywords: Topological insulator; proximized superconductivity; Green functions; odd-frequency topological 
regimes; bound states; conductance spectra 


2 



The search for a qubit that is robust against local decoherence has led to extensive 
studies of the so-called Majorana zero-energy mode in materials with strong spin-orbit cou¬ 
pling (SOC) combined with s-wave superconductors. Nanowires and topological insulators 
are promising candidate materials with strong spin-orbit coupling [1-15]. Experimentalists 
have already reported signatures of the Majorana zero-energy mode [16-18], where zero-bias 
conductance peaks are the main features observed in this context. 


Besides the possibility to host Majorana modes, topological materials are also interesting 
for the study of unconventional p-wave superconductivity by itself, p-wave superconductivity 
gained renewed interest after the prediction that the p-wave pairing symmetry in He^ would 
lead to half-quantum vortices with potential application to the held of quantum computing 
[19-21]. Sr 2 Ru 04 is believed to have p-wave symmetry and together with the prediction of 
the existence of a nodal gap it has led to an extensive study of this material [22-29]. Also 
superconductor/ferromagnet/superconductor junctions are studied for their prospective to 
switch from a dominant s-wave state to a p-wave state [30-35]. Also these devices have 
potential application to quantum computation, through which the understanding of the 
behaviour of p-wave superconductivity has become a held on its own. 


In this paper we study the s and p-wave correlations that exist in three-dimensional topo¬ 
logical insulators (3D TIs). We hrst discuss in which regime a dominant p-wave correlation 
is present and under which conditions a dominant p-wave correlation can lead to a zero- 
energy Majorana bound state. We derived the expressions for the bound states formed at 
a 3D TI in the “Majorana” regime. We show the relation between odd-frequency pairing 
and Majorana zero-energy modes by using a Green function approach. We use those new 
insights to calculate the conductance spectra of proximity induced superconducting 3D TIs 
with and without broken time-reversal symmetry. The main focus is to understand how the 
modelled tunneling conductance spectra arise from the s and p-wave correlations. We will 
see that the combination of a zero-bias conductance peak together with conductance dips at 
the gap energy are distinguished features for p-wave correlations. We demonstrate that it is 
not necessary to be in the “Majorana” regime in order to observe signatures of the p-wave 
correlations. 
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I. PAIRING WAVE FUNCTION AND MAJORANA-MODES IN A 3D TI 


In order to determine the pairing relations for a material with strong SOC we start with 
the relation [E — H]G = 1 where G is the Green function of the system, 1 the identity matrix 


and H is in the basis 
given by 


V’l, V’j) V’f) W [6]- For a system with spin orbit coupling [E — El] is 


E -\- (^p ~\~ M —— iXyky 

- \xkx iXyky E €p — M 


0 A 

-A 0 


0 

A 


E — Cp — M 

~Xxkx tXyhy 


-A 

0 

Xxkx tXyky 

E — ep + M 


( 1 ) 


where ep = fi— -—. Here, Xx,y is the SOC strength in x and ^/-direction. At the end we will 
set the SOC strength equal to the Fermi velocity and ep = fi to discuss the superconducting 
correlations in the topological surface states, fi is the chemical potential, A is the induced 
superconducting gap and M is the energy of a perpendicular magnetic held, either externally 
applied (Zeeman term) or due to an exchange interaction. Note, that the Zeeman contri¬ 
bution of a magnetic held can be considerable, due to the large g^-factor of most relevant 
materials. We neglect the orbital contribution of the applied held. 

By taking the inverse of this matrix equation we can obtain the Green’s function G, 
expressed as 


Gii Gi2 
G 21 G 22 

The elements Gij are blocks of 2x2 matrices. The diagonal blocks describe the propagation 
of the electrons and holes separately. The oh-diagonal blocks describe the combined electron 
and hole propagation, i.e. they describe the propagation of the Cooper pairs in the s-wave 
proximized topological insulator which is the quantity of interest in this paper. When we 
represent 


the Green function is given by 


G{k,E) = 


E-H{k) 



{a — h d ^ c) ^ {c — db ^ d) ^ 
{b — d c~^ d)~^ {d — c d~^ b)~^ 


( 3 ) 

( 4 ) 
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We find 


Gi 2 —i— [{E'^ — + ‘^{—^xkx^p + iXykyM)ai 

-2{-Xykyep - iX^k^M)a2 - 2MEaz] a2, 

Z = — (E‘^ — ~ 4:{~Xxkx€p + iXykyM)"^ 

— A{—Xykyep — iXxkxMY — 4:E‘^M‘^, 


Bk,M —€p + + {XxkxY + i^ykyY ~ ( 7 ) 

To discuss the symmetry of the anamolous Green functions, we use the Matsubara rep¬ 
resentation which can be obtained by analytical continuation, E + i6 ^ iun, 

Gi2 ~ + ‘^{—Xxkx^p + iXykyM)ai 

-2{-Xykyep - iXxkxM)a2 - i2Mund-3] (J 2 ( 8 ) 

=iifo + f ■ cr)d 2 , 

Z = — {—cOn ~ Bk^nvi)"^ ~ ^{~Xxkx^p + iXykyM)'^ 

- 4{-Xykyep - iXxkxMf + (9) 

fo{k, iu) = - Bk,M), (10) 

fi{k,iu) = 2^{-Xxkxep + iXykyM), (11) 

A 

f2(^kjiu!^ ^y^y^p 'iXxkxJXI^^ ( 12 ) 

fz{k,iuj) =-2i^Mujn- (13) 

The total wave function of a pair of fermions should be asymmetric and the total wave 
functions can be described by the product of an orbital (or parity), spin and frequency 
term. Even-frequency pairing means that a function is even in Un- If we consider singlet 
(which is an odd function under spin permutation) s-wave (orbitally symmetric) pairing, 
the pairing wave function E{= G 12 ) should satisfy the relation Es{k) = Es{—k) in order for 
the wave function to be antisymmetric when the pairing is even in frequency. For p-wave 
triplet pairing the relation Ep{k) = —Ep{—k) applies [36]. 
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We note that Z is an even function of a;„ and is an even-parity function. The singlet 
component /o belongs to the even-frequency even-parity class (ESE). The two equal-spin 
components (/i and ) are even-frequency spin-triplet odd-parity class (ETO). Finally, /s 
belongs to odd-frequency spin-triplet even-parity class (OTE). 

For a 3D topological insulator and the SOC strength in x and y-direction in Eq. (1) 
are equal to fi and the Fermi velocity v respectively. The matrix [E — H] in the relation 
[E — H]G = 1 is then given by 

E + jj, — M —v\k\E^ 0 —A 

—E + a + M A 0 

. , (14) 

0 A E — fj, + M —v\k\e 

—A 0 —v\k\e^^ E — fi — M 

where M is a magnetization term and 6 is the angle between kx and ky. It is instructive to 
compare the 3D TI with semiconductors with strong spin-orbit coupling. Although there 
are analogies between them, also several differences exist between them. For comparison, 
we show therefore the model for nanowires in the supplementary material. 

We assume in this section that fi > Edp where DP stands for the Dirac point, but the 
opposite “hole” regime of fi < Epp can also easily be obtained. For a chemical potential 
far above the Dirac point we can make additionally the assumption that /i 3> A. If we also 



Fii ~ 2A\k\vE^ (M + fi) /Ztd 
Fij ~ A (M^ - \kf - fi^) /Zti, 

Fji ~-A {M^-\kfv^-f,^) /Zpi, 

Fjj ~ 2A\k\ve-^^ (M - fi) /Zti, (16) 

Zpi = —B^j — 4 {—vkxfi — ivkyM)^ — 4 {—vkyfi + ivkxM)^ , 

BIj = {fA + v‘^\k\^ -M^Y + 2A2 {fA + v^\k\^ - M^) . (17) 

G 21 is related to G 12 by complex conjugation. The factor shows the chiral p-wave 
character of an s-wave proximized 3D topological insulator. 
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For the moment we assume that we have even-frequency symmetry, [f^ in Eq. (15) can 
be neglected in the regime of low energy excitations. The appearance of non-neglible odd- 
frequency symmetry terms will be discussed in section III.) The second and third relation 
in Eq. (16), therefore, correspond to s-wave pairing and the other two with p-wave pairing 
relations. 

The energy dispersion relation of the s-wave proximized topological insulator can be 
obtained by diagonalizing the corresponding Hamiltonian as described above, and is found 
to be 

E = ±\JEf ± 2a/E^ (18) 

where Ef = and Eg' = In the limit of 

/i ^ A this can in good approximation be written as 

E = ±/i ± ^/v^\k\^ + M\ 

Now we have all the formal work done to consider the pairing symmetry in three different 
regimes: \M\ = 0, M < /i and the “Majorana” regime. 

A. Regime \M\ = 0 

Assume a homogenous magnetization term (for example caused by an external magnetic 
held). In the regime \M\ /i we have the situation as depicted in Fig. 1. There is a gap 
opening indicated by Ai at positive momentum. At negative momentum a gap with the 
same size is opening. Using the approached dispersion relation we have n = v\k\ at the gap 
opening. Substituting this relation into Eq. (16) we obtain the following anamolous Green 
function relations at Ai 

E^^r^2A\k\\kf\EyZTi, 

_P|j ~ —2A\k\‘^ /Zti, 

Fjl ~ 2A\k\‘^ /Zti, 

Eii^-2A\k\\kf\e-^^/ZTi. 

We see from the above that the matrix elements corresponding to the s-wave pairing are 
in magnitude as large as the matrix elements corresponing to the p-wave relations. So in 
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Momentum 


FIG. 1. Energy dispersion relation for a 3D topological insulator in the regime \M\ <C /r. The solid 
line corresponds to \/Et' — 2-^Eg/ and the dashed line to \/ Et' + 2yJEsi. The graph is plotted for 
A in the order of /r for clarity. 


the case of neglible M we have an equal admixture of s-wave and p-wave correlations. This 
result is similar to the result found by Tkachov in Ref. [37, 38] for a topological insulator. 
In Ref. [15] Fu and Kane transform to another basis, Ck = + in which the 

Hamiltonian becomes equivalent to a spinless + Wy (dominant p-wave). Although this 
transformation is an elegant way to show that p-wave relations are present, when doing 
predictions for the experimental outcome in this basis, a tranformation also has to be made 
to the s-wave superconductor deposited on top. This would also cause an additional phase 
shift in the normal superconductor through which one can then conclude that p-wave is not 
longer dominant. Nonetheless, the existing proposals of Fu and Kane for Majorana devices 
[15] still hold as it was already noted that breaking time-reversal symmetry can make Px+ipy 
dominant, as we will see next. 





B. Regime M < n 


As we turn on the perpendicular magnetization term, it follows from the dispersion 
relations that we have at Ai, . The anamolous Green function relations 

become 

Fii ~ 2Ae*"-^(/i - M) (/i + Mf/Zn, 

Tu ~ —2A 

Tji ~ 2A /Ztii 

Fjj ~ -2Ae-'®\/(/i - Mf (/i + M)/Zti. 

As soon as time-reversal symmetry is broken the increase in M causes the p-wave compo¬ 
nent to become larger in magnitude than the s-wave component. This can be intuitively 
understood by notifying that a magnetic held aligns the spins parallel to the held. Therefore 
spin triplet pairing is favoured above spin singlet pairing. It is therefore possible to already 
observe dominant p-wave features in the regime M < p. 


C. The “Majorana” regime 

Unfortunately, the presence of dominant p-wave correlations still does not mean that a 
Majorana zero-energy mode exists as soon as we apply magnetization, because the zero- 
energy mode is not yet fully localized. As already note by Fu and Kane [15], the 2D TI 
surface has no edge, which prevents a localized Majorana bound state to form. The same is 
true in a 2D topological insulator where ID edge states can go around the 2D topological 
insulator. We therefore need a diherent device in order to localize the zero-energy mode. 

One way to create an edge and break time-reversal symmetry is shown in Fig. 2(a). The 
M denotes the magnetization induced by the ferromagnet. We assume here that we have 
an insulating ferromagnet so that the current is only going through the surface states. A 
different way is to have a Josephson junction with a phase difference between the two super¬ 
conductors [15]. For now, we focus on the device in Fig. 2(a). In order to have a localized 
zero-energy mode, the Majorana mode needs to be localized both at the superconducting 
side and at the ferromagnet side. At the superconducting side the superconducting gap 
is doing the job. At the ferromagnet side, the magnetization has to be large enough such 
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(a) 


(c) 



FIG. 2. (a) A ferromagnet insulator(FI) can both break time-reversal symmetry and localize 

the Majorana zero-energy mode in a TI. The location of the Majorana mode is indicated by 
the white arrow, (b) The dispersion relation of the topological insulator at the ferromagnet side 
for a magnetization smaller and larger than the chemical potential. In the former, the chemical 
potential still lays in the Dirac cone. In the latter, the states around the chemical potential are 
gapped out and the zero-energy mode is localized at the interface between the superconductor and 
ferromagnet. (c) Calculation of the surface Andreev bound state. To model the suppression of 
the superconducting gap of the superconducting topological insulator (STI) we assume a normal 
part in the order of the coherence lenght at the interface between the superconductor and the 
ferromagnet. 


that the Fermi level is inside the gap (see Fig. 2(b)). Then, the zero-energy mode is fnlly 
localized. To be more specihc, in order to have a proper localized Majorana mode one has 
to satisfy the relation M{x) > + jj? at the ferromagnet side. We also inclnded the 

snpercondncting gap as a fnnction of position to take the proximity effect into acconnt. 
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II. SURFACE ANDREEV BOUND STATES AND MAJORANA ZERO-ENERGY 


MODES 


So far we have looked at properties of the bulk, i.e. the inner part of the s-wave proximized 
2D surface states of the topological insulator. We are now going to discuss the properties at 
the interface between the superconductor and ferromagnet where the Majorana zero mode is 
formed in the Majorana regime. It is expected that different symmetries are present at the 
surface or edge of a superconducting system as noticed in Ref. [39, 40]. Due to a hnite size, 
translational symmetry is broken in that direction. The orbital symmetry (even or odd) is 
then not a well-dehned parameter anymore, giving rise to mixed odd and even symmetry at 
the interface. The spin pairing symmetry, however, is not affected by the hnite size. The 
only option, then, to have an asymmetric wave function is to allow for odd-frequency states 
at the interface. 

Figure 2(a) shows the superconductor/ferromagnet device needed to localize a Majorana 
state. In experiments it is useful to use a ferromagnetic insulator (FI) to ensure that the 
magnet itself does not shunt the device [41]. We will, therefore, indicate the ferromagnet by 
FI in the remaining part of the paper. 

Assume, that we are in the regime with a fully localized Majorana state, i.e. M > /i. 
The resulting SABSs can be calculated by considering a small non-superconducting region 
in the order of the coherence length to model the suppression of the gap near the edge (Fig. 
2(c)). We use the wave functions described in our previous work (Ref. 42) to calculate the 
SABSs. The direction of the transmitted waves in the superconductor is chosen such that k is 
conserved parallel to the interface. We assume that the electrons in the non-superconducting 
part are conhned and therefore, that they completely reflect back at a; = —L. 

The phase difference that is picked up in one round trip (Fig. 2(c)) is 


so 


2L{\ke\ — \kh\) — 29 + 371 + 2 arccos E/A = 27rn 
ALE 

= - -- — 29 + Sti + 2 arccos E/ A 

nvf cos 9 


(19) 


E 

A 


cos 


^—37r /2 + 9 + 
= — sin (9 + 


2LE \ 
hvfcos 9 J 
2LE \ 
hvf cos 9 J 


( 20 ) 

(21) 
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Taking the limit L —)■ 0 we have 


( 22 ) 

(23) 


^ - siii(S) 

When the magnetization of the FI is in the exact opposite direction, the resulting ABS is 
^ = sin( 6 *). This is exactly the SABS of a chiral p-wave superconductor [43]. 

In order to determine the features of the SABS and its relation to Majorana zero en¬ 
ergy modes, we consider the regimes around E = Q and i? 7 ^ 0 in more detail in the 
following. 


A. Around E = 0 


In Fig. 2(c) we can calculate the reflection and transmission coefficients of the electron 
and holes at the interface by using continuation of the wave function. From these coefficients 
it follows that the total wave function in the superconducting part in the electron-hole basis 
V’fct > V’H, V’l fct ’ V’l is given by 

f (E- tn) \ 
e-iO(E - in) 

A 

( (E + zQ) \ 

{E + in)E^ 

—Ae*® 

A 


ie 




^-ikr^-rl^O 


! 


E + iVt 


\ 


^ikr^-rlio 


(24) 


/ 


where hi = — E"^ and .^0 = VA"^ — E'^/vp- Around E = 0 (and so around 0 = 0) the 

expression for the Majorana zero mode becomes 

/ cos (/ex -I-7r/4) ^ 

2gW4 gQg _ 7 ^/ 4 ^ 

2gW4 gQg -f 7r/4) 

y 26 “*^^/^ cos {kx — 7 r/ 4 ) J 


T = g-^Ao 


(25) 
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Note, that this expression is consistent with the result found by Tanaka and Asano in 
a semiconductor/superconductor device [6]. In their case, an additional magnetic held is 
applied through which only the spin up part survives of Eq. (25). 

To give a complete picture of the properties of these Majorana zero-energy mode, we cal¬ 
culated the frequency symmetry of this zero-energy state. The anamolous Green’s function 
can be obtained by the relation 


F{E,r,r'-,e)u 


E + i6 — En E + i6 + En 


(26) 


where u and v are the electron and hole parts respectively. 
The relationship at x = x' = 0 satishes 


F/3,„(r, F, p) = - [F„,^(r, -F, p)]* (27) 

where F denotes the conjugated of F and a and /3 are the spin-indices[36], which means the 
Cooper pairs are odd-frequency pairs. The spatial symmetry of Eq. (25) satishes the s-wave 
symmetry. Therefore the boundary of a pure p-wave even-frequency superconductor is a 
pure s-wave, odd-frequency state. The appearance of odd-frequency states near zero-energy 
is also noted in calculations on and px + ipy superconductors. [44] 


B. F / 0 

In a similar way as discussed for zero energy we can calculate the anomalous Green 
function for energies larger than zero. We hnd that there is an admixture between even and 
odd-frequency. The odd and even-frequency parts are, respectively, given by 


F„,, ^ Vl - FVA2, 

(28) 

Eeven ^ F/A. 

(29) 


Near zero-energy, the odd-frequency part is dominant and the even-frequency contribution 
grows with F. The emergence of a (partly) odd-frequency amplitude at the surface/interface 
due to spatial non-uniformity [39, 45-47] or interorbital pairing [48] is well-known for su¬ 
perconducting systems. However, in this case, we have a full odd-frequency state at zero 
energy. Therefore the conclusion of Asano and Tanaka [6] holds that in these devices pure 
odd-frequency Gooper pairs and the Majorana zero-energy mode are one and the same thing. 
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(a) 


(b) 



FIG. 3. (a) In the calculation of the conductance we consider two electrodes: the proximized TI 
by the FI and the proximized TI by the s-wave superconductor, (b) In order for the Andreev 
formalism to hold, the distance between the FI and the superconductor should be smaller than the 
mean free path and the coherence length. 

III. TI/STI TUNNELING CONDUCTANCE 

In the study of the conductance we consider the system where one electrode consists 
of the topological insulator proximized by the FI and the other electrode of the topological 
insulator proximized by the s-wave superconductor as depicted in Fig. 3(a). We assume here 
that the magnetization is smaller than the chemical potential in the magnetic topological 
insulator electrode. We modelled this magnetic TI/STI interface by using the wave functions 
in the magnetic TI and STI as described in Ref. [42] and match them at the interface such 
that 


i’in + aiJh + blpe = tei^Se + thi’Sh 


(30) 


where i/jh is the reflected hole, -i/e the reflected electron, i/jse the transmitted quasi-particle 
in the electron branch of the superconductor, 'ipsh the transmitted quasi-particle in the hole 
branch of the superconductor and a, 6, te and th are the Andreev, normal, electron trans¬ 
mission and hole transmission coefficients respectively. We calculated then the conductance 
by using the relation 


7r/2 

/ d9cos{9)T{e,eV) 
G —7r/2 


Go -/2 

/ d9cos{9)T{9,eV)\A=o 

-7r/2 


(31) 


where T {9,E) = 1 + |ap — and 9 is the angle between kx and ky. The Andreev formalism 
holds when the distance between the superconducting and the ferromagnet side is smaller 
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than the coherence length and the mean free path {Imean) as drawn in Fig. 3(b). We choose 
here to divide the conductance by the energy dependent normal conductance to ensure that 
the normalized conductance goes to one at large voltage. 

We hrst consider the case with time-reversal symmetry, i.e. without the FI on top of the 
TI. By means of a gate electric held, the chemical potential at the superconductor side can 
be tuned independently from the non-superconductor side. Therefore, we make from now on 
a distinction between the chemical potential at the superconductor and non-superconductor 
side. In Fig. 4(a) we plotted the normalized conductance for different ratios of the STI 
chemical potential, /is, and TI chemical potential, jiTi- A conductance peak at el^ = A is 
appearing for larger barriers due to Andreev resonance at this energy, similar to the s-wave 
case. Opposite to the s-wave case, the value of the conductance never drops to zero, even for 
a large mismatch where the barrier height goes practically to inhnity = 1/100). The 

conductance inside the gap is close to one instead which can be understood by considering 
different angles of incidence. The large barrier causes the propagation direction in the 
superconductor to be almost equal to 0 = 0 as follows from the conservation of parallel 
momentum: ds = arcsin [sin 6 ht i/f^s)- The spin in the topological insulator with 6 = 0 
sees no barrier because of Klein tunneling (which gives perfect Andreev reflection for all 
energies) and for 6 = j2 the spin mismatch is the largest which results in the lowest 
transparency (perfect normal reflection). For angles inbetween, the transparency of the 
barrier increases continuously from one to zero for increasing angle of incidence. The angle 
averaged conductance is therefore close to one. 

The presence of a non-zero conductance for energies below the gap, even for large barriers, 
distinguishes the proximized TI surface from conventional s-wave symmetry superconduc¬ 
tors. Although the p-wave pairing is not dominant in the time symmetric situation, the 
p-wave correlations are encoded in the non-zero conductance for energies below the gap in 
the presence of a barrier. The fact that the gap does not go to zero is due to the spin- 
momentum locking which is also responsible for the p-wave correlations. 


IV. TI/STI TUNNELING CONDUCTANCE WITH BROKEN TIME-SYMMETRY 

To modulate the TI/STI tunneling conductance with broken time-reversal symmetry we 
take the set-up as displayed in Fig. 3(a). We will use a relatively large chemical potential 
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(a) 


(b) 





FIG. 4. (a) Conductance spectra for a TI/STI configuration for different mismatches in chemical 
potential, (b) Conductance spectra for a magnetization M/A=95, ^'r7/A=100 and different mis¬ 
matches in chemical potential, (c) Conductance spectra for a magnetization M/A=10, /rT//^=100 
and different mismatches in chemical potential, (d) The conductance spectra for different values of 
magnetization and fiTi/^ = /.is/A=100. Even if the magnetization is smaller than the chemical 
potential, clear signatures of the SABSs are still visible. 


compared to the superconducting gap to model the experimentally realistic regime. The 
topological insulators so far have the Dirac point close or at the same energy as the bulk 
bands so as long as the chemical potential is in the gap, the chemical potential is likely much 
larger than the superconducting gap. We argued in earlier work[42], that a magnetization 
in the case of a ferromagnet insulator will only be a few percent of the chemical potential. 
However, to see the effect of the broken time-symmetry we also modulated the conductance 
spectrum for a magnetization value of 0.95fiTi- This result is shown in Fig. 4(b). There are 
two distinguished features compared to an N/S interface, namely a zero-bias peak (ZBP) and 


16 






















a conductance dip at eV = A. The appearance of these two features for large magnetization 
is also reported by Linder et a/. [41], Both features are due to the formation of SABSs at the 
interface. In fact, the SABSs have a constant density of states at energies between —A and 
+A as follows from Eq. (22), but from conservation of parallel momentum it follows that the 
angles near zero have a larger transparency. The coherence peaks at the gap energy in the 
low barrier limit, are replaced by resistance peaks at the same energy. The coherence dips 
at A arise from the spectral weight transfer from high to low energy due to the formation 
of SABSs. This results in a decrease of the density of states at A. Therefore, next to a zero 
bias conductance peak, also coherence dips are a distinguished feature of a p-wave order 
parameter compared to s-wave. 


Note, that the surface states are not fully gapped yet with a magnetization of 95% of the 
chemical potential, i.e. the chemical potential is not inside the gap at the FI side. We are 
therefore not in the topologically non-trivial regime. However, the increased reflecion at the 
interface or enhanced p-wave pairing symmetry due to magnetization, causes the SABS to 
form at the interface. If we use a magnetization value of 10% of the chemical potential in 
the TI, Fig. 4(c) is obtained. We see there is only a small effect on the conductance spectra 
around zero energy at A compared to the time-reversal symmetry case in Fig. 4(a). Figure 
4(d) shows the effect of increasing the magnetization for fiTi = hs' = lOOA. Gradually, 
the conductance at zero energy increases while at A the conductance decreases due to the 
formation of SABSs for larger magnetization. It shows that the ZBP and conductance dip at 
the gap energy come together. That is, if one has a dominant p-wave pairing wave function 
and if one observes a ZBP, then also the conductance dips are present. As soon as the ZBP 
becomes less pronounced, as is the case for lower magnetizations, the conductance dip at 
eV = A is less pronounced. The observation of both a ZBP and a conductance dip is a 
strong signature for a p-wave state. For practical ferromagnets that have a magnetization in 
the order of a few percent of the chemical potential of the TI, it is not expected to see this 
particular p-wave signatures back in the conductance spectra. A pronounced deviation from 
s-wave behaviour starts to occur when M/A = 60 (see Fig. 4(d)). With a typical exchange 
energy of 40 meV [49, 50] it means that one has to gate tune the chemical potential to about 
70 meV from the Dirac point, which seems doable experimentally. [51, 52] 
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V. DISCUSSION AND CONCLUSION 


An equal admixture of s and j9-wave correlations exists in a 3D topological insulator 
proximized by an s-wave superconductor. By inducing a perpendicular magnetization, the 
p-wave pairing becomes dominant. For a magnetization energy as large as the chemical 
potential, the topological regime is entered with Majorana modes. 

We studied the symmetry of the bound state of a 3D topological insulator in the Majorana 
regime. Green function techniques show us that the pure odd-frequency pairing state at 
interfaces is equivalent to a zero-energy Majorana mode. 

In the conductance spectra the increasing p-wave pairing can be observed by the presence 
of a ZBP together with a conductance dip at the gap energy. In the time-reversal symmetric 
case, the non-zero sub-gap conductance in the presence of large interface barriers, indicates 
that p-wave correlations are present. 
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the Dutch Foundation for Fundamental Research on Matter (FOM), the European Research 
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